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REFERENCE:  Contract  NObs.  86437 
SF- 013- 11-01, 

Task  1359 


TO:  Chief,  Bureau  of  Ships 

Department  of  the  Navy 
Washington  25,  D.  C. 

ATTN:  Code  345 

FROM:  Edwin  J.  Waller,  Project  Director 

School  of  Civil  Engineering 
Oklahoma  State  University 
Stillwater,  Oklahoma 

SUBJECT:  Final  Report 


This  report  summarizes  the  work  performed  under  contract 
NObs  86437  during  the  period  February  1,  1962, to  January  31,  1963. 
Work  is  continuing  under  contract  NObs  882  97  which  terminates  Jan  ¬ 
uary  31,  1964. 

The  over-all  objective  of  this  project  was  to  aid  the  Bureau 

of  Ships  in  prevention  of  excessive  liquidborne  noise  levels  at  critical 

) 

points  in  shipboard  piping  systems.  This  entails  the  following: 

A.  To  develop  a  standard  method  of  evaluating  the  noise 
producing  characteristics  of  pumps  and  flow  control  devices 
in  piping  systems. 

B.  To  develop  the  necessary  information  leading  to  the  pre¬ 
paration  of  standards  and  specifications  for  pumps  and 
control  valves  by  studying  the  response  of  various  typical 
piping  arrangements  to  noise  inputs  typical  of  the  pumps 
and  control  devices  and  presenting  the  results  in  terms 
of  the  acoustical  levels  at  the  discharge  of  the  piping  sys¬ 
tem. 

C.  To  carry  out  necessary  mathematical  analyses  as  a  part 
of  B.  above  which  yield  pertinent  information  useful  for 
improving  design  and  analysis  techniques  for  the  over-all 
pump-piping  system  including  design  of  noise  filters. 

Liquidborne  noise  is  generated  in  pumps  and  control  valves  by 
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fluid  turbulence,  cavitation,  and  by  the  acceleration  of  the  liquid  by 
moving  mechanical  parts  such  as  vanes  or  pistons.  The  noise  thus 
created  at  the  pump  or  valve  is  transmitted  (i.e.  propagated)  in  the 
form  of  pressure  pulsations  throughout  the  piping  system.  Additional 
noise  is  generated  by  the  fluid  as  it  flows  through  the  piping  system, 
depending  upon  its  configuration.  This  project  was  concerned  only 
with  the  noise  which  is  generated  at  the  pump  or  control  valve. 

The  noise  intensity  at  any  given  point  in  the  system  depends 
on  the  dynamic  response  capability  of  the  piping  system  and  the  inten¬ 
sity  of  the  noise  generated  by  the  pump  or  control  device.  There  are 
two  fundamental  means  of  altering  the  noise  level.  These  are  either 
to  redesign  the  piping  system, minimizing  the  noise,  or  to  decrease  the 
noise  level  at  the  source. 

Generally  the  piping  system  and  the  pumps  are  designed  by 
different  persons  for  performance  criteria  which  gives  little  or  no 
consideration  of  noise  generation  and  subsequent  transmission.  In 
some  systems  the  noise  level  is  acceptable  while  in  others  it  is  not. 
This  is  due  to  differences  in  the  piping  system  parameters  such  as 
elasticity,  viscosity,  and  density  of  the  liquid;  the  elasticity,  diameter, 
and  length  of  the  pipe;  and  the  terminating  conditions  of  the  piping 
system.  Thus,  the  parameters  of  the  piping  system  have  considerable 
influence  on  the  noise  spectrum. 

Figure  la  shows  the  types  of  discontinuities  encountered  in 
shipboard  piping  systems.  These  include  changes  in  size  (A),  looping 
of  lines  (B  to  C),  and  stub  lines  (D).  Additional  complications  arise 
when  multiple  sources  and  multiple  terminal  points  are  looped  into 
one  system,  as  indicated  by  Figure  lb. 


4 


Specific  Objectives  and  Resume  of  Results. 

A  literature  survey  of  pertinent  publications  and  reports  was 
the  first  objective  of  this  project.  Material  was  supplied  by  the  Okla¬ 
homa  State  University  Library  and  from  the  Bureau  of  Ships.  The 
results  are  shown  in  the  Bibliography  of  this  report..  Although  the 
results  of  the  survey  were  informative,  no  practical  method  of  deter¬ 
mining  the  response  of  complex  piping  systems  was  found  in  the  litera¬ 
ture. 

The  next  objective  was  to  establish  a  practical  method  of  de¬ 
termining  the  response  of  the  complex  piping  system  by  mathematical 
analysis.  This  analysis  was  carried  out  and  presented  with  The  Third 
Quarterly  Report  of  this  project.  Necessary  revisions  have  been  made 
and  are  presented  in  Chapters  I  to  IV  of  this  report.  The  results 
include  the  mathematical  methods  used  in  the  analysis,  the  idealiza¬ 
tions  involved  in  the  physical  model  of  the  piping  system,  the  techni¬ 
ques  and  formulas  for  finding  the  transfer  functions  relating  the  pres¬ 
sure  spectra  at  various  points  in  the  system,  and  general  problems 
involving  input  specification.  The  results  are  formidable,  but  indicate 
a  technique  of  analysis  which  can  be  carried  out.  The  limitation,  is 
computer  facilities. 

Considerable  time  was  spent  in  developing  and  evaluating 
various  analysis  methods  with  the  objective  of  finding  the  most  feasible 
procedure.  The  most  promising  of  these  is  discussed  in  Chapter  V. 
Although  the  transfer  functions  developed  in  Chapter  III  are  valid  for 
all  frequencies,  a  substantial  amount  of  information  could  be  gained 
from  the  maximum  and  minimum  of  the  transfer  functions  if  they  could 
be  determined.  However,  the  results  shown  in  Chapters  V  and  VI 
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giving  upper  and  lower  bounds  of  pressure  spectra  are  unreasonably 
extreme  and  impractical  for  use  in  a  simplified  analysis. 

Other  results  are  shown  in  Chapter  VI.  These  include  transfer 
function  calculations  for  suction  lines  in  a  system  similar  to  the  auxi¬ 
liary  sea  water  system  of  the  USS  TINOSA  (  SS(N)606)  in  which,  allow¬ 
able  pressure  spectra  are  found  for  the  frequency  range  0  <  f  <  10,  000 
cps.  The  results  "de -bugged"  part  of  a  general  program  which  will  be 
used  in  future  studies.  Essentially  this  technique  involved  calculation 
of  the  transfer  functions  in  complex  number  form. 

The  mathematical  technique  under  development  at  Oklahoma 
State  University  is  being  adapted  for  computer  solution.  The  programs 
developed  will  use  realistic  piping  and  flow  data  obtained  from  the  auxi¬ 
liary  sea  water  system  of  the  TINOSA.  During  the  last  two  months  of 
the  contract  period,  some  time  .was  spent  in  collecting  information  such 
as  attenuation,  characteristic  impedance,  and  phase  constants  for  the 
sea  water  system  from  the  drawings  of  the  system,  and  in  analysis  of 
heat  exchangers  to  determine  their  dynamic  response.  Work  continues 
toward  the  over -all  goals. 


CHAPTER  I 


TRANSFORM  METHODS  IN  SYSTEM  ANALYSIS 

1.  1  Introduction. 

Operational  calculus,  or  transform  methods,  has  proven  to  be 
an  elegant  tool  in  the  analysis  of  linear  systems.  The  system  under 
consideration  is  no  exception.  Certain  relationships  involving  trans¬ 
form  methods  are  dealt  with  in  this  chapter. 

1.2  Fourier  Series. 

Adequate  description  of  the  dependent  variable(s)  involved  in  an 
analysis  is  important.  For  periodic  functions,  this  is  often  accom¬ 
plished  by  the  Fourier  series  or,  as  to  be  used  here,  the  complex 
form  of  the  Fourier  series. 

Although  any  independent  variable  can  be  used,  the  dependent 
variable  will  be  considered  a  function  of  time  f(t)  and  periodic  in  time 
T.  Then*108*,* 

00  •  x 

f(t)  =  -T^F(Wn)e  n  *  j  =a/71~  f1-1) 

n=-oo 

where  the  complex  coefficient 

/T/2  .  . 

-jw  t 

f(t)  e  n  dt  ,  (1.2) 

-T/2 

and  the  nth  harmonic  of  the  fundamental  angular  frequency 

*  Numbers  in  parentheses  refer  to  references  in  Bibliography. 
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(1.3) 


1.  3  Fourier  and  Laplace  Transforms. 

Mathematical  definitions  and  justifications  of  Fourier  and  Laplace 
transforms  are  available  in  many  sources  The  purpose  of 

this  section  is  to  establish  the  interrelationship  of  the  Fourier  series 
description  of  a  function  and  the  Fourier  transform. of  the  function 
Equation  (1.1)  can  be  written 


1  'V 

fw  =  TT  1  P(“n>  e  “1 


(1.4) 


n=-oo 


where  w1  is  the  fundamental  angular  frequency.  As  the  period  T  tends 
to  infinity,  the  function  f(t)  tends  to  be  aperiodic,  the  fundamental  an¬ 
gular  frequency  becomes  a  differential  of  angular  frequency  dto,  and 
the  nth  harmonic  angular  frequency  con  becomes  a  continuous  angular 
frequency  to.  In  the  limit,  f(t)  becomes  aperiodic,  and  Equation  (1.4) 
becomes 


f(t) 


1 

2  7T 
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F(w)  e^wt  dw 


(1.5) 


where 


F(«) 


/ 


f(t)  e  ^  dt 


(1.6) 


from  Equation  (1.‘2) 

Equations  (1.5)  and  (1.6)  are  recognized  as  a  Fourier  transform 
pair.  The  corresponding  Laplace  transform  relationships  are 


e+joo 

/ 


e-jco 


f(t) 


1 

2?rj 


F(s) eSt  ds, 


(1.7) 


8 


and 

oo 

A 

F(s)  = 

j  f(t)  e"st  dt 

(1.8) 

W 

0 

where  s  is  the  complex  variable  a  +  jo. 

If  the  function  f(t)'is  defined 

as 

f(t)  =  <j 

[f(t>,  t>0 
[_  0  ,  t  <  0  , 

(1.9) 

it  is  evident  that  the  Laplace  and  Fourier  transform  pairs  becpme  iden¬ 
tical  when  s  is  replaced  by  jo  . 

/I  nn\ 

1.4  Spectrum  Analysis  K±'J  . 

From  the  results  of  preceding  discussion,  it  is  concluded  that  a 
dependent  variable  can  be  described  in  either  the  time  or  frequency 
domain.  In  the  time  domain,  the  variable  can  be  expressed  by  an  ap¬ 
propriate  analytical  function  or  by  a  Fourier  series  if  periodic.. 

The  frequency  domain  description  of  a  periodic  function  is  evident 
from  Equation  (1.  1).  The  function 


T 


F(o, 


n 


for  n  =  0,  +1,2,3...  gives  a  measure  of  the  periodic  component  of 
f(t)  at  the  nth  harmonic  frequency  wn>  In  general,  F(on)  is  complex, 
and 

1  111  i®n 

T  FK]  =  T  |  F(wn}  |e  (L10) 

where-  ©n  is  the  phase  (angle  of  F(on).  Plotting 


1_ 

T 
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versus  w  (or  n)  results  in  a  discrete  spectrum  known  as  the  amplitude 
spectrum  of  f(t).  Similarly,  the  graph  6 ^  versus  u>n  (or  n)  is  the  phase 
spectrum. 

For  the  non-periodic  function  f(t),  the  function 

-5J-  F(»)  -  ~  I  F(,|  t  eie<“>  (1.11) 

is  a  continuous  function  of  w  and  is  in  general  complex.  It  is  the  com¬ 
plex  continuous  spectrum  of  f(t).  The  amplitude . density  spectrum  of 
f(t)  is 


and  its  phase  density  spectrum  is  0(w).  Since  division  by  2ir  in  Equation 
(1.  11)  is  simply  a  normalizing  factor,  the  function  |F(w)|  shall  here¬ 
after  be  referred  to  as  the  amplitude  density  spectrum  of  f(t). 

1.  5  Transfer  Functions  and  Frequency  Response. 

The  analysis  of  a  linear  system  by  Laplace  transform  methods 

(99  114) 

often  results  in  a  relationship  of  the  form  '  ’ 

F0(s)  =  H(s)F.(s)  (1.12) 

where  Fq(s)  and  F.(s)  are  the  Laplace  transforms  of  the  system  output 

fQ(t)  and  input  L(t),  respectively.  The  function  H(s)  is  usually  referred 

to  as  a  transfer  function  or  system  function. 

If  f.(t)  is  defined  as  in  Equation  (1.  9),  the  variable  s  can  be  re- 

* 

placed  by  jw  in  Equation  (1.  12)  to  give 


*When  s  is  replaced  by  jw  in  a  function  F(s),  the  result  will  be 
denoted  F(w)  rather  than  F(jw). 
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FoM  =  H(u)F.(u).  (1.13) 

The  function  Fq(io)  will  be  called  the  frequency  response  of  the  system. 

The  amplitude  density  spectrum  of  the  system  output  in  terms 
of  the  input  spectrum  and  system  function  follows  by  taking  the  absolute 
value  of  both  sides  of  Equation  (1.  13).  Thus, 

|FqH|  =  [H(<o)|  |Fi(co)f  .  (1.14) 

The.  term  transfer  function  will  be  used  interchangeably  to  des¬ 
cribe  either  H(s),  H(oj),  or  |H(co>|  . 


CHAPTER  II 


PHYSICAL  SYSTEM 


2,  1  Introduction. 

This  chapter  contains  the  following: 

1.  Derivation  of  differential  equations  describing  compressible, 
turbulent  flow  of  water  or  similar  liquid  in  a  non-rigid  cylin¬ 
drical  conduit. 

2.  Discussion  of  conditions  imposed  on  system. 

3.  Laplace  transform  domain  solution  of  describing  equations. 

4.  System  frequency  response  and  discussion  of  system  para¬ 
meters. 

The  results  of  this  investigation  comprise  a  major  part  of  the 
necessary  information  required  to  determine  the  response  of  a  piping 
system  such  as  the  one  shown  in  Figure  la. 

2.2  Derivation  of  Describing  Equations. 

Figure  2  depicts  a  liquid-filled  cylinder  with  non-rigid  walls 
under  compressible,  turbulent  flow  conditions.  Conditions  imposed 
on  the  system  follow  and  hereafter  will  be  referred  to  by  number: 

1.  The  flow  is  one -dimensional. 

2.  The  frictional  resistance  is  a  function  of  some  nonlinear 
operation  on  the  fluid  velocity  and  system  parameters.  It 
can  be  expressed  as  the  product  of  a  constant  and  the  instan¬ 
taneous  flow  rate  variation. 
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3.  All  energy  dissipation  is  accounted  for  by  Condition  2. 

4.  The  rate  of  change  in  fluid  velocity  with  respect  to  length 
is  negligible  as  compared  to  its  rate  of  change  with  respect 
to  time. 

5.  The  instantaneous  pressure  is  constant  over  the  cross 
section. 

6.  The  rate  of  change  in  mass  density  with  respect  to  length  is 
negligible  as  compared  to  its  rate  of  change  with  respect  to 
time. 

7.  The  velocity  of  wave  propagation  in  the  liquid  cylinder  is 
constant. 


\ 

Figure  2 . 


The  total  velocity  u^.  of  a  liquid  particle  is  considered  to  be  the 
sum  of  a  mean  velocity  u  and  an  instantaneous  velocity  variation  u 
assumed  positive  in  the  direction  of  mean  flow,  i.e. ,  u^  =  u  +  u  . 
Likewise,  the  total  pressure  p^.  is  the  sum  of  a  mean  pressure  p  and 
an  instantaneous  pressure  variation  p,  i.  e. ,  p^  =  p  +  p. 

Considering  conditions  1,  2,  3,  4,  and  6  and  assuming  the  positive 
lateral  direction  x  to  oppose  the  mean  flow,  the  following  relationships 
will  completely  describe  the  system 
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Equation  of  motion: 


&u,  ,  9p, 

—  -  -  —■  +  Guf  =  0 

9t  p  9x  t 


(2.  1) 


where  p  is  the  liquid  mass  density  and  G  is  a  constant. 
Equation  of  continuity: 

•3u, 


-p 

9  9x  9t 


0  . 


(2.2) 


Equation  of  state: 


9p, 


K' 


9  P 


(2.3) 


where  K'  denotes  the  bulk  modulus  of  the  liquid  and  conduit  combined. 


It  is  given  by 


(114) 


K'  = 


KbE 
KD  +bE 


(2.4) 


where  K  denotes  liquid  bulk  modulus,  b  is  conduit  wall  thickness,  E 
is  the  modulus  of  elasticity  of  wall  material,  and  D  denotes  inside  dia¬ 
meter. 

Equations  (2.  2)  and  (2.  3)  are  combined  to  give 


^t 

8x 


K’ 


9pt 

9t 


-  0  . 


(2.  5) 


Integrating  Equations  (2.  1)  and  (2.  5)  over  a  control  volume  A6x  and 
using  condition  5, 


and 


4* 


+ 


0  , 


aqt  .  A  9Pt 
dx  K'  9t 


0 


(2.6) 


(2.7) 
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where  A  denotes  pipe  area  and  B  is  a  constant.  The  total  volume  flow 
rate  is  the  sum  of  a  mean  flow  rate  q  and  an  instantaneous  flow  rate 
variation  q  assumed  positive  in  the  direction  of  mean  flows  i.  e. , 

=  q  +  q.  The  nonlinear  equation  (2.6)  is  linearized  as  follows^). 

If  the  flow  rate  is  considered  constant  at  a  cross -section)  Equation 
(2.  6)  becomes  the  slope  of  the  pressure  grade  line,  i.  e. , 


•o  — : n 
B  q 


fpq  2 

2  DA2 


(2.8) 


where  f  is  the  Darcy-Weisbach  friction  factor.  Then 


B 


fpq  2;n 

2  DA2 


(2.9) 


where  pf  denotes  the  pressure  needed  to  overcome  the  frictional  re¬ 
sistance  of  the  pipe.  The  value  of  n  is  estimated  from  experience  or 
experiment  and  should  range  between  1.65  and  2.  05, 

The  nonlinear  term  in  Equation  (2.6)  is  rewritten  and  expanded 
as  follows: 


Bq“  =  B'(  q  +  q)n 

=  Bqn(l  +  i)n 

-^[l+n(3)+i!|pil  (  3  . ; 


(2.10) 


The  series  in  Equation  (2.  10)  is  convergent  for  'q  <  q  ,  and 
since  q«q  ,  the  series  is  assumed  sufficiently  approximated  by  its 
first  two  terms.  Equation  (2.6)  is  now 

np. 


8x 


+ 


Ill  + 

6t  i 


+ 


iq 


q  =  0 


(2. 11) 


15 


The  desired  describing  differential  equations  for  instantaneous 
variations  in  flow  rate  and  pressure  follow  directly  from  Equations 
(2.  11)  and  (2.  7)  by  replacing  q^  and  p^  by  their  respective  sums.  They 
are: 


Ox 


JL  il 
a'  at 


4- 


npf 

~?W 


q 


0  ; 


(2. 12) 


and 


dq  ,  A  8p 

aF  +1T  at 


o  . 


(2.13) 


2.  3  Imposed  Conditions. 

The  conditions  imposed  on  the  physical  system  are  discussed  and 
validated  to  varying  degrees  in  the  following. 

1.  In  any  real  piping  system,  a  propagated  wave  is  never  one¬ 
dimensional,  because  the  elasticity  of  the  pipe  walls  allows 
radial  motion.  But  over  a  certain  frequency  range  from  0 
to  radians  per  second  there  exists  only  one  propagated 
wave,  or  one  mode  of  propagation,  which  is  the  axial,  or 
„  plane  wave,  mode.  This  mode  is  denoted  the  (0,  0)  mode. 

The  cut-off  frequency  wc  at  which  the  first  radial  mode  of 
propagation  occurs  j^the  (0,  1)  mode[]  for  nondissipative  con¬ 
ditions  can  be  determined  by  the  following  equations  given 
by  Jacobi 


JG(z)  2p2b 

zJ^z)  P^D 


2Eb 

■p.  2  2  * 

P1Dc  z 


(2.  14) 


and 


w 


c 


2c 

D 


z 


01 


(2.  15) 


i 
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■where  z  is  some  characteristic  value,  JQ  and  J  are  Bessel 
functions  of  the  first  kind,  and  and  p ^  are  the  mass  den¬ 
sities  of  the  fluid  and  pipe  material,  respectively.  The  first 
real  characteristic  value  that  satisfies  Equation  (2.  14) 
is  substituted  into  Equation  (2.  15)  to  obtain  w  ,  For  8"  stan¬ 
dard  steel  pipe,  f  =  6900  cycles  per  second.  For  smaller 
sizes  of  pipe  the  value  of  coc  increases  substantially.  To  the 
best  knowledge  of  the  authors,  the  subject  of  cut-off  fre¬ 
quencies  for  dissipative  conditions  is  nonexistent  in  the  litera¬ 
ture.  The  cut-off  frequency  should  be  somewhat  higher  for 
dissipative  conditions  than  that  given  by  Equations  (2.  14)  and 
(2.  15).  It  is  concluded  that  for  a  relatively  rigid  pipe,  the 
one -dimensional  condition  is  valid  for  frequencies  up  to  the 
cut-off  frequency. 

2.  The  perturbation  process  used  in  the  derivation  of  the  des¬ 

cribing  equation  to  linearize  the  frictional  dissipation  term 
was  proposed  and  verified  experimentally  for  low  frequencies 
by  Waller  There  seems  to  be  no  good  reason  for  dis-  ' 

counting  its  validity  in  the  frequency  range  of  one -dimen¬ 
sional  propagation. 

3.  The  effect  of  heat  transfer  on  a  propagated  wave  could  be 

considered  resulting  in  a  third  describing  differential  equa- 

(15) 

tion  as  done  by  Brown  v  ,  but  it  is  believed  that  the  rather 
empirical  frictional  dissipation  term  used  in  the  describing 
equations  will  account  for  all  energy  dissipation. 

4.  The  condition  that  the  rate  of  change  in  fluid  velocity  with 
respect  to  length  is  negligible  has  been  used  for  many  years 
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in  describing  wave  propagation  in  liquid  cylinders.  Phillips^6^, 
in  showing  the  acoustic  wave  equation  adequate  for  describing 
turbulent  flow  in  liquid-filled  steel  pipes, has  shown  this  term 
to  be  negligible. 

5.  If  the  one -dimensional  wave  condition  is  valid,  it  follows  that 
the  instantaneous  pressure  is  essentially  uniform  over  the 
cross  section. 

6.  As  in  Condition  4,  the  rate  of  change  in  mass  density  with 
respect  to  length  has  classically  been  considered  negligible. 

7.  Strictly  speaking,,  the  velocity  of  wave  propagation  in  a 
liquid  cylinder  with  non-rigid  walls  is  not  constant  with  fre¬ 
quency.  With  increasing  frequency,  the  flexural  vibrations 

of  the  walls  tend  to  impede  the  propagation  causing  a  decrease 
in  the  wave  velocity.  Mathematical  expressions  to  determine 
the  velocity  of  wave  propagation  as  a  function  of  frequency 
are  given  in  the  paper  by  Jacobi^41\  For  a  relatively  rigid 
pipe,  such  as  a  standard  steel  pipe,  the  velocity  of  wave 
propagation  in  the  pipe,  a,  is  essentially  constant,  and^^ 


a 


(2.  16) 


Rather  than  calculating  the  velocity  of  propagation  by  Equa¬ 
tions  (2. 4)  and  (2.  16),  it  can  be  determined  for  different 
pipe  sizes  and  fluids  from  nomographs  available  in  many 
papers  and  books 


2.4  Solution  of  Differential  Equations. 

The  discussion  is  facilitated  by  introducing  three  system  -para¬ 


meters.  They  are: 


Coefficient  of  inertia, 
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Coefficient  of  resistance, 

r 

1  q 

Coefficient  of  capacitance, 


C 


(2.17) 


(2.  18) 


(2.19) 


In  terms  of  these  parameters,  the  describing  differential  Equa¬ 
tions  (2.  12)  and  (2.  13)  are: 

-  -8-g -fe-  ?)--  +  L  *)•■  +  R  q(x,  t)  =  0  ;  (2.  20) 

and 

-  dS-feAi  +  C  -9P- .fet)  =  o  .  (2, 21) 

3x  dt  ' 

Using  Laplace  transform  methods,  Equations  (2.20)  and  (2.21) 
are  solved  for  pressure  and  flow  rate  in  the  transformed  time  domain 
as  follows 

Transforming  Equations  (2.  20)  and  (2.  21)  with  the  Laplace  inte¬ 
gral 

f(t)  e"st  dt  , 

+  (R  +  sL),  Q(x,  s)  p  0  ,  (2.22) 

+  s  CP(x,  s)  =  0 


00 

F(s)  = 


9P  (x,  s) 
9x 


and 

9Q  (x,  s) 
ax 


(2.  23) 
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where  p(x,  0)  and  q(x,  0)  are  zero.  Transforming  the  above  equations 
with  thb  Laplace  integral 


*<*)  -  f 

J  n 


f(x)  e  ^  dx  , 


and 


-  XP(X,  s)  +  (R  +  sL)  Q  (X,  s)  =  -P(0,s),  (2.24) 

-  X^(X,s)  +  sC  £  (X,  s)  =  -  Q(0,  s)  .  (2.25) 

Solving  Equations  (2.  24)  and  (2.  25)  for  ^(X,  s)  and  4(X,  s), 


and 


where 


and 


P(X)S)  =  ZeQ(0.s)  (  g  y-2  )  +P(0,s)(— n  -X--  2'  V  (2.  26) 
VX-y  /  V  X  -  y  / 


§(x,s)  =  Q(°.*i(^'.X  a) 12,2,1 


y  =  sC  (R  +  sL)  , 


,2  R  +  sL 
Jc  sC 


(2.28) 


(2.  29) 


Evaluating  Equations  (2.  26)  and  (2.  27)  with  the  inversion  integral 


f(x)  = 


.e=joo 

1  1  Xt 


27rj 


/ 

e-joo 


£(X)  e''-1’  dX 


and 


P(x,  s)  =P(0,  s)  cosh  yx  +  Z£  Q(0,  s)  sinh  yx  ,  .  (2,30) 


Q(x,  s)  =  Q(0,s)  cosh  yx  +_-pl^  s). 

c 


sinh  yx  . 


(2.31) 
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Equations  (2.  30)  and  (2.  31)  provide  the  Laplace  transform  domain 
solution  for  pressure  and  flow  rate  fluctuation  at  a  point  in  a  pipe  in  terms 
of  a  propagation  coefficient  y,  a  characteristic  impedance  Z  ,  and  pres¬ 
sure  and  flow  rate  variation  at  a  reference  point  x  =  0. 

If  the  mean  flow  is  in  the  'direction  of  positive  lateral  distance, 
it  is  an  easy  matter  to  verify  that  the  preceding  analysis  is  valid  if  x  is 
replaced  by  -x  in  all  functional  entities.  Thus,  if  mean  flow  is  from  ' 
right  to  left  in  Figure  2, 


P(x,  s)  =  P  (  0,  s)  cosh  yx-Z^  Q(0,  s)  sinh  yx  ,.  (2,  32) 

and 

Q  (,x,  s)  =  Q  (  0,  s)  cosh  yx  -  s)  sinh  yx  .,  (2.33) 

c 

This  will  continue  t<5.  be  true  in  all  further  developments. 

The  foregoing  analysis  is  based  on  turbulent  conditions,  but  it 
holds  equally  well  for  laminar  flow.  For  laminar  flow. 


32p  iq 
AD2 


and  n  =  1.  Thus, 

32  ju 


R 


AD 


(2.34) 


(2.  35) 


where  p  denotes  dynamic  viscosity.. 

For  a  pipe  having  one  end  closed  the  mean  flow  is  zero.  Here 
the  analysis  breaks  down  only  in  the  determination  of  the  frictional  re¬ 
sistance  term.  It  appears  reasonable  to  assume  the  coefficient  of  re¬ 
sistance  for  this  case  to  be  the  same  as  for  laminar  flow. 
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2.  5  Frequency  Response  and  System  Parameters. 

The  variables  p(x,t)  and  q(x,t)  are  considered  the  output  of  the 


system.  The  input  is  p(0, t)  and  q(0,  t).  Then,  the  frequency  response 
of  the  system  is  obtained  from  Equations  (2.  30)  and  (2.  31)  by  replacing 


s  by  jw. 

That  is, 

P(x,w)  =  P(0,w)  cosh  yx 

+  ZcQ(0,w)  sinh  yx. 

(2. 36) 

and 

Q(x,w)  =  Q(0,w)  cosh  yx 

+~(  ^ sinh  yx 

(2. 37) 

where 

and 

7 2  =  jw  C  (R  +  jwL)  , 

c 

(2.  38) 

2%  _  R  +  ioj  L 
c  jwC 


(2.  39) 


The  propagation  coefficient,  y,  can  be  expressed  in  the  form 
a  +  jj3  by  finding  the  roots  of  Equation  (2.  38).  The  attenuation  constant, 
a  ,  is  given  by 

l 

a  =  ~  [Jr2  +  w2L2  -  0)Lj  J 

The  phase  constant,  j3,  is  given  by 

o  =  [if  . 

The  characteristic  impedance,  Z ^  ,can  also  be  expressed  in 
complex  form  by  finding  the  roots  of  Equation  (2.  39).  It  is  given  by 


(2.40) 


(2.41) 


(0  "  j'a)  • 


(2.42) 


22 


Expanding  Equation  (2.  40)  in  powers  of 


1 

2 


(2.43) 


If  the  above  series  is  assumed  sufficiently  approximated  by  its  first 
term,  Equation  (2.  43)  becomes 


arv 


R 

2La  ’ 


(2.44) 


By  a  similar  argument.  Equation  (2.  41)  becomes 


j3~ 


co 

a 


(2.45) 


Assuming  nto  be  small  as  compared  to  jS,  the  characteristic  impedance 
can  be  considered  real  and  given  as 

Zc~La.  (2.46) 

Actual  calculation  of  a,  j8,  and  Z£  for  various  systems  by  both 
exact  and  approximate  formulation  has  shown  the  approximations  to  be 
of  sufficient  accuracy  with  the  possible  exception  of  very  low  frequen¬ 
cies  Therefore,  Equations  (2.44),  (2.45),  and  (2.  46)  will  be  used 

in  further  developments. 


CHAPTER  III 


TRANSFER  FUNCTION  DETERMINATION 

3'.  1  Introduction. 

Using  Equations  (2.  36)  and  (2.  37),  transfer  functions  relating 
output  pressure  and  flow  to  input  pressure  are  obtained  for  a  single 
pipe  and  parallel  pipes.  Extension  of  the  single  pipe  analysis  to  pipes 
in  series  is  described.  The  transfer  relations  for  single  and  series 
pipes  are  obtained  in  two  different  forms,  both  giving  identical  results 
and  each  having  merit  for  different  situations.  The  development  is 
based  on  total  length  of  a  component,  but  can  easily  be  generalized 
for  any  arbitrary  point  in  the  system. 

3.  2  Single  Pipe. 

Figure  3  depicts  a  single  pipe  which  is  the  ith  component  of  a 
system.  The  receiving  and  sending  ends  of  the  pipe  are  denoted  by  i 
and  i+1,  respectively.  The  complex  pressure  and  flow  rate  variations 
at  a  point  i  are  denoted  P^  and  Q„  The  absolute  values  of  these  vari¬ 
ables  are  and  <3^  . 

The  impedance  Z.  at  a  point  i  is  defined  as  the  ratio  of  pressure 
to  volume  flow  rate  at  that  point,  i.e. , 


P. 
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The  absolute  value  of  Z.  is  Z^,  and  its  phase  angle  is  0.,  i.  e. , 

_  30i 

Z.  =  Z.e  1  .  (3.2) 


P.., , 
i+l 


Q* 


i+l 


x. 

i 


P. 

i 

Qi 


flow 


#i 


i. 

1 


i+l 


Figure  3.  Single  Pipe. 


Denoting  y  and  Z^  of  pipe  i  as  y ..  and  Z^,  Equations  (2.  36)  and 
(2.  37)  can  be  manipulated  to  yield 


and 


P.x1  =  P 

i+l 


.[ 


cosh  +-£ 


Z  . 
ci 


sinh  y  .1. 

'  i  l 


=  P.G.  , 

i  i 


('3.  3) 


S+l  =  P 


*[+ 


cosh  y .£.  + 

1  T  1 


i  i  Z  . 

Cl 


sinh  y 


A] 


=  P.H. 

l  l 


(3.4) 


where  and  IT  are  transfer  functions.  The  complex  functions  G^ 
and  EL  can  be  expressed  in  more  suitable  form  for  numerical  use  by 
replacing  by  a.  +  j/L  .  Equations  (3.  3)  and  (3.  4)  represent  one  form 
of  transfer  function  determination. 

The  reflection  coefficient  at  a  point  i  is  defined  as 


Z.  -  Z  . 

4  1  Cl 

z.  +  z  . 

l  ci 


(3.5) 
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The  magnitude  of  is  P  ,  and  its  phase  is  0. .  Using  Equation 
(3.5),.  Equations  (2.  36)  and  (2.  37)  are  manipulated  to  give 

-2y,.£. 


P..,  =  P.e 
l+l  i 


Yiii  r  i  +  1 1 1 1 


-n  i' 


(3.6) 


ano- 


y.J.  f\  T\  2Vi 


i  i  1  -  Fe 

o  =  p  5 _ :  j 

Wi+1  i  Z  . 

1  + 


Cl 


(3.7) 


The  second  useful  form  of  transfer  function  determination  follows 
from  taking  the  absolute  value  of  both  sides  of  Equations  (3.  6)  and  (3.  7) 
This  gives 

p  =  R  [1  +  Ai  +  2Ai  cos^.'j  \ 


i+1 


1  +  p2  +  2  P  cos  ej 

*  1  4  1  1 


-  pi  Gi  • 


(3.8) 


and 


a.l .  f1jA2  ..  ")  1 

=  p  i  i  1  +  A.  -  2A  cos  ip, 

l+l  l 


'i  2 


Z  .  ]  *-»2  p 

C1  1  1  +  I  .  +  2  I  .  cos  0. 

1  l  *  l  l 


=  P.  H. 

1  l 


(3.  9) 


where 


r-2a.£. 

,e  1  1 

L 


(3.  10) 


and 


i p.  =  -2BJ.  +  0.  . 
1  11  1 


(3.  11) 
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The  amplitude  density  spectrums  of  the  pressure  and  flow  rate 
at  the  point  of  i  +  1  are  obtained  either  by  taking  the  absolute  value  of 
results  from  Equations  (3.  3)  and  (3.  4)  or  by  employing  Equations  (3.  8) 
and  (3.  9). 

• 

3.3  Series  Pipes. 

It  has  been  shown  that  transfer  functions  can  be  obtained  relating 
the  output  at  a  point  in  a  series  system  to  the  input  at  any  other 
point  The  resulting  functions  are  quite  complicated  and  seem 

to  possess  no  advantages  over  the  following  metho'd  of  analysis. 

Consider  the  series  piping  system  shown  in  Figure  4.  The  no¬ 
tation  follows  that  of  the  preceding  section.  The  logical  method  of 
analysis  follows  from  either  Equations  (3.  3)  and  (3.  4)  or  Equations 
(3.8)  and  (3.9). 


flow, 


n+1 


Figure  4.  Series  Piping  System. 


The  analysis  starts  at  the  point  1  where  end  conditions  are  known, 
i.  e.  ,  P  ^  and  are  known.  Using  Equations  (3.  3)  and  (3.  4)  or  Equa¬ 
tions  (3.  8)  and  (3.  9),  P^  and  or  P2  and  are  obtained.  The  pro¬ 
cedure  is  repeated  for  each  segment  until  the  desired  point  is  reached. 


27 


To  repeat  the  above  procedure,  the  impedance  Z..  of  each  ith 
segment  must  be  determined.  If  Equations  (3.  3)  and  (3.  4)  are  being 
employed, 

Vrr  ■  <3- 12> 

1-  1 

With  reference  to  Equations  (3,  8)  and  (3.  9)  , 

Z.  =  ~-1—-  e^1  (3.  13) 

1  Hi-1 

where 

2A.  .  sin  il/.  , 

0.  =  tan'1  - 1-7 1  -g-  (3.14) 

1  1 -Af 1 

i-l 

as  can  be  shown  by  manipulating  Equations  (3.  6)  and  (3.  7). 

3.  4  Parallel  Pipes. 

A  parallel  piping  system  is  shown  in  Figure  5.  The  parallel 
components  are  not  necessarily  dimensionally  or  materially  identical. 


m* 

# 

m 

m 

1 _ 

1 

..  _ _ _ J 

— 

- 1 

1 - ■/  Ui — — - ' 

1  A  A  A  / - - -1 

_ 

II. 

~{A  — 

2.' 

#2 

2 

flow 

■i+1 

■ 

i 

1' 

*1 

1 

Figure  5.  Parallel  Piping  System 


The  development  starts  with  P^  and  Z^  assumed  known  which  in 
turn  determines  Q.^.  The  following  relationships  are  used: 
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p.  = 

1 

P1  = 

P2  =  •  •  • 

=  Pm  ’ 

(3.  15) 

Qi  = 

Q  + 

q2  +  .  .  . 

+  Qm  j 

(3.  16) 

Pi+1 

=  pi< 

=  p2,=  ■ 

‘  *  =  Pm'  ; 

(3.  17) 

Qi+1 

=  Ql' 

+  q2,  +. 

•  •  +  Qm'  ; 

(3.  18) 

Using  the  notation 

Bk  =  cosh  Yk  ik  ,  (3.19) 

Tk  =  sinh  Ykik  ,  (3.20) 

and 

Mk  =  sinh  y  i  (3.21) 

ck 

where  k  refers  to  pipe  number.  Equations  (2.  36)  and  (2.  37)  give 

pk-  =  pkBk  +  QkV  (3-22) 

and 

Qk'  =  <\Bk  +  PkMk  *  <3-23> 

k  =  1,2,...,  m.  Using  relationships  (3.  15)  and  (3.  17),  subtraction 
of  Equation  (3.  22)  for  k  =  1  from  each  of  the  remaining  Equations  (3.  22) 
yields  the  m  - 1  equations 

TlQl-TkQk  =  Pi(Bk-Bl>  *  <3-24) 

k  =  2,3,...,  m.  Adding  Equations  (3.23)  and  substituting  relation¬ 
ships  (3.  15)  and  (3.  18)  gives 

m  m 

«i+i  1  BA  -  pi  2  Mk  ■  <3-  25> 

k=l  k= 1 

Equations  (3.  16),  (3.24),  and  (3.  25)  provide  m+1  linearly  inde¬ 
pendent  equations  in  the  m+1  unknowns  and  Qk>  k=l,  2, . . .  ,  m. 
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Thus,  Q^+1  can  be  determined,  and  once  any  is  found.  Equation 
(3.22)  yields  P^+^  directly.  Also,  the  pressure  and  flow  rate  at  any 
point  in  a  pipe  k  can  be  determined  by  finding  and  using  Equations 
(2.  36)  and  (2.  37) . 

The  results  of  the  analysis  can  be  expressed  in  the  general  forms 


P.  .  -  =  P.  G. 

l+l  i  i 


and 


(3. 26) 


Vi  =  Pi  Hi 


(3. 27) 


where  G.  and  H.  are  transfer  functions.  For  m=2, 

ii' 


T1  [fB2-Bl>Zi+TJ 


Gi  B1  +  Z.  (T1  +  T2) 


(3.  28) 


and 


H. 


.  t  <B2  ‘  El>2  +  <“l  +  M2><V  T2)]  Z1+  BiT2+  B2Tj 


Zi  <V  T2> 


(3.29) 


The  functions  G^  and  H.  become  increasingly  complicated  as  m  in¬ 
creases. 

If  the  parallel  pipes  are  dimensionally  and  materially  alike. 


*  Q,  =—  Q.  , 
k  m  i 


(3.  30) 


and 


Qk’  Qi+1 


(3.31) 


Then,  Equations  (3.22)  and  (3.23)  give 

f  1  Tk 
P..  «  =  P-  [  B,  +  — —  TB— - 
l+l  i  \  k  m  Z. 


=  P.  G.  , 
i  i 


(3.  32) 
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and 

■  pi  ( %  +  ”Mk) 

(3.  33) 

=  P.  H. 
i  i 

for  any  k  =  1,2, ... ,  m. 

3.  5  Summary. 

The  results  of  this  section  make  it  possible  to  determine  pres¬ 
sure  and  flow  rate  fluctuation  at  any  point  of  a  system,  such  as  shown 
in  Figure  la,  in  terms  of  an  input  pressure  and  impedance  at  some  other 
point  of  the  system.  The  determination  of  this  input  pressure  and  im¬ 
pedance  is  discussed  in  the  next  chapter. 

The  transfer  functions  and  usually  the  input  involved  in  the  analy¬ 
sis  of  a  system  are  frequency  dependent.  The  output  of  the  system  is 
then  dependent  on  frequency;  in  fact,  it  varies  considerably  over  re¬ 
latively  small  frequency  increments.  Thus,  analysis  of  a  piping  system, 
for  a  large  range  of  frequencies  involves  an  impractical  amount  of  nu¬ 
merical  computation.  For  this  reason,  the  possibility  of  finding  the 
maximum  and  minimum  system  output  is  investigated  in  Chapter  V. 


CHAPTER  IV 


INPUT  DETERMINATION 

4.  1  General. 

Application  of  analysis  methods  developed  thus  far  is  dependent 
on  ability  to  determine  an  input  pressure  and  impedance  at  some  point  ' 
in  a  piping  system.  In  the  general  sense,  the  input  to  a  system  is 
caused  by  a  system  element  producing  a  flow  disturbance.  Typical 
flow  disturbing  elements  are  pumps  and  valves. 

If  the  input  is  deterministic.,  i.e.  ,  expressible  as  a  function  of 
time,  the  Fourier  transforms  of  functional  relationships  for  pressure 
and  flow  rate  fluctuation  can  usually  be  obtained  This  yields  the 

necessary  information  concerning  input  for  analysis.  Two  determinis¬ 
tic  inputs  often  encountered  in  piping  systems  are  reciprocating  pumps 
(94,  95)  and  valve  clogures  (76,  93). 

In  the  case  of  centrifugal  or  similar  type  pumps,  the  input  is 

non -deterministic,  or  random.  Many  other  flow  disturbing  devices 

such  as  valves  produce  random  inputs.  Statistical  methods  must  then” 

(20  36  92) 

be  introduced  in  an  analysis  incorporating  such  inputs  '  '  ’  ' . 

Statistical  analysis  lends  itself  toward  uncertainty,  because  determin¬ 
ation  of  statistical  properties  concerning  pump  and  valve  inputs  is  at 
best  questionable. 

The  foregoing  discussion  of  input  determination  has  been  directed 
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toward  a  problem  of  prediction,  i.  e.  ,  given  a  system  containing  a 
flow  disturbing  element,  what  is  the  pressure  and  flow  rate  variation 
at  some  other  point  in  the  system  ?  Another  problem  of  equal  impor¬ 
tance  is  one  of  specification,  i.e. ,  given  a  system,  what  is  the  allow¬ 
able  pressure  and.  flow  rate  variation  at  a  point  in  the  system  so  as 
not  to  exceed  some  specified  variation  at  another  point?  More  speci¬ 
fically,  given  a  system,  what  is  the  allowable  pressure  amplitude  den¬ 
sity  spectrum  at  a  point  in  the  system  so  as  not  to  exceed  some  speci¬ 
fied  pressure  amplitude  spectrum  at  a  termination  point  of  the  system 

The  methods  of  analysis  developed  in  Chapter  III  are  easily  applied 
to  the  problem  of  specification  for  systems  with  deterministic  or  random 
inputs.  The  remainder  of  this  chapter  and  Chapter  V  is  directed  toward 
problems  of  this  type. 

Assuming  a  pressure  amplitude  density  spectrum  is  specified  at  a 
termination  point  of  a  system,  the  impedance  at  that  point  must  be  de¬ 
termined  to  employ  the  transfer  relations  developed  in  Chapter  III.  This 
impedance  is  denoted  by  Z^.  Closed  and  open  end  terminations  are  dis¬ 
cussed  in  the  following., 

4.  2  Closed  End  Impedance. 

For  a  pipe  with  one  end  rigidly  closed,  the  flow  rate  fluctuation  at 
the  closed  end  is  zero.  The  closed  end  impedance  is  then  infinite.  It 
is  also  easily  shown  that  the  reflection  coefficient  is  1. 

A  stub  line  such  as  shown  in  Figure  la  as  pipe  D  is  a  case  of 
closed  end  termination.  In  this  case,  the  impedance  of  the  stub  line 
at  the  junction  is  necessary  in  analysis  of  the  complete  system.  From 
Equations  (c>.‘3)  and  (3.4),  the  impedance  is  calculated  directly  by 
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Z 


=  Z 


cD  coth'7DiD 


(4.  1) 


4.  3  Open  End  Impedance. 

The  impedance  at  an  open  end  termination  of  a  piping  system  is 
usually  considered  zero  on  the  hypothesis  that  pressure  fluctuation  at 
that  point  is  zero.  This  assumption  is  valid  only  for  very  low  frequen¬ 
cies  where  the  range  of  validity  depends  on  pipe  size. 

If  the  velocity  at  the  open  end  is  considered  harmonic  and  of  con¬ 
stant  profile,  the  problem  becomes  the  classical  acoustics  problem  of 
determining  the  impedance  of  a  rigid  circular  piston  set  in  an  infinite 
baffle.  The  foregoing  assumptions  seem  reasonable  and,  at  least, 
offer  a  simple  method  of  determining  the  open  end  impedance. 

An  excellent  development  of  the  solution  to  the  afore  mentioned 
problem  can  be  found  in  the  text  by  Kinsler  and  Frey  The  de¬ 

sired  result  is 


PC 

A 


•[r^x) 


+  3  Xx(x) 


(4.2) 


where  c  is  the  velocity  of  wave  propagation  in  the  'infinite  medium,  and 


x 


ojD 

c 


(4.3) 


The  area  of  the  open  end  is  A,  and  the  diameter  is  D.  Also, 


Rx(x)  =  1 


2Jx(x) 


x 


2>4  2.  42  •  6 


2-  42- 62- 8 


(4.  4) 


and 


2K  (x) 
Xl(x)  =— i— 

X 
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where  K^x)  is  a  related  Bessel  function.  For  small  values  of  x., 

2 

Rl(x)-~  ,  (4.6) 

and 

Xl»~fr  .  (4.7) 

The  functions  Rj(x)  and  X^(x)  are  shown  graphically  and  tabulated 
numerically  in  Kinsler  and  Frey's  book.  The  graphical  representation 
is  reproduced  in  Chapter  VI,  Figure  8.  The  functions  are  also  approxi¬ 
mated  for  computational  purposes  in  that  chapter. 

Although  the  open  end  impedance  depends  on  frequency,  it  can  be 
considered  essentially  constant  over  frequency  bands.  This  fact  is  used 
in  the  next  chapter. 


CHAPTER,  V 


MAXIMUM  AND  MINIMUM 
ALLOWABLE  PRESSURE  SPECTRUM 

5 .  1  Intro  duction . 

An  allowable  pressure  level  at  some  point  in  a  piping  system  for 
a  given  frequency  can  now  be  determined.  To  obtain  a  complete  allow¬ 
able  pressure  spectrum  for  a  sizable  frequency  range,  a  large  number 
of  calculations  must  be  made.  For  this  reason,  the  problem  of  deter¬ 
mining  a  maximum  and  minimum  allowable  pressure  spectrum  is  now 
approached. 

5.2  General. 

Systems  to  be  studied  will  have  open  end  termination  as  discussed 
in  Section  4.  3  of  this  report.  This  type  of  termination  is  normally  the 
case  encountered  in  sea  connected  systems.  The  specified  pressure 
amplitude  spectrum  at  the  termination  is  assumed  constant  over  a  fre¬ 
quency  band.  This  will  allow  easy  adjustment  to  the  true  overboard 
specifications  of  pressure.  Thus,  the  specified  pressure  spectrum  and 
impedance  at  the  termination  are  considered  constant  over  a  frequency 
band  where  the  bands  cover  the  frequency  range  of  interest,  e.  g.  ,  one- 
third  octave  bands  from  0  to  10,  000  cps. 

The  object  of  the  investigation  is  then  to  determine  the  maximum 
and  minimum  allowable  pressure  levels  at  some  point  in  the  system 
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for  each  frequency  band.  Once  these  levels  are  determined,  the 
following  remarks  can  be  made  about  the  actual  pressure  amplitude 
occurring  at  the  point  for  a  given  frequency  in  a  band.. 

If  the  pressure  amplitude  is  greater  than  the  maximum  level, 
the  pressure  amplitude  at  the  termination  is  greater  than  the  specified 
amplitude.  If  the  pressure  is  less  than  the  minimum  level,  the  termi¬ 
nation  amplitude  is  less  than  the  specified  pressure.  If  the  pressure 
is  between  the  maximum  anc}:  minimum  levels,  the  termination  pres¬ 
sure  can  be  either  above  or  below  the  specified  pressure  amplitude. 
For  this  case,  the  allowable  pressure  for  the  frequency  in  question 
must  be  calculated  to  reach  a  conclusion. 

If  these  maximum  and:  minimum  levels  can  be  determined,  the 
problem  of  specification  is  reduced  considerably.  Much  effort  by  the 
authors  has  been  devoted  to  determining  these  extremes.  It  has  been 
concluded  that  the  maximum  and  minimum  levels  cannot  be  determined 
analytically  for  any  but  the  simplest  system,  i.  e. ,  a  single  pipe. 

The  next  section  describes  a  method  for  analytically  determining 
the  desired  extremes  for  a  single  pipe.  Section  5.  4  discusses  the 
problems  encountered  in  attempts  to  extend  the  analysis  to  pipes  in 
series.  The  developments  are  based  on  simple  arguments.  Methods 
of  variational  calculus  attempted  provided  no  information  other  than 
the  obvious. 


Equations  (3.8) 

(5.  1) 


5.3  Single  Pipe. 

The  single  pipe  system  is  shown  in  Figure  6. 
and  (3.  9)  give 


P2  =  Pre 


ffjiljl  +Ai  +  2A^  cos 


1 
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and 
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p  _L_ 
r  Z 


l  1  Jl  +  A2  -  2A^  cos  4/ 


cl  1  + 


+  2 


cos  6 


(5.2) 
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Figure  6.  Single  Pipe. 


Since  all  end  conditions  are  assumed  constant  over  a  frequency 
band,  the  only  frequency  dependent  parameter  in  Equations  (5.  1)  and 
(5.2)  is  /3.p  or  ip  ^  =  -  20^  +  @r  .  Obviously,  the  maximum  P2  occurs 
when  cos  ^  =  1,  and  the  minimum  P2  occurs  when  cos  i jj  ^  =  -1.  Simi¬ 
larly,  cos  4/ ^  =  1  gives  minimum  Q2,  and  cos  41  ^  =  -1  gives  maximum 


<^2  • 


The  desired  maximum  and  minimum  pressure  levels  for  a  fre¬ 


quency  band  are  then 


Pr  eVl  (1+A1) 


2max 


rr 


T'*2  +  2  r*  cos  0 
r  lr 


(5.3) 


and 


Pp  eal£l  (1-AX) 


'  2  min 


”l  +  F2  +  2F  cos  el-i 
_  1  r  *r  j 


(5.4) 


The  number  of  calculations  required  to  cover  the  frequency  range  of 
interest  depends  only  on  the  number  of  bands  chosen  to  cover  the  range. 
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For  a  narrow  bandwidth,  there  is  a  possibility  the  values  of 
cos  Tp-y  used  to  obtain  Equations  (5.  3)  and  (5.  4)  do  not  occur  in  the  band. 
If  this  discrepancy  cannot  be  allowed,  it  is  an  easy  matter  to  determine 
the  maximum  and:  minimum  values  of  cos  P  ^  for  the  band  and  use  them, 
in  Equations  (5.  1)  and  (5.2). 

5,4  Series  Pipes. 

Two  pipes  in  series  are  shown  in  Figure  7.  The  maximum  and 
minimum  allowable  pressure  levels  for  a  particular  frequency  band  at 
point  2  are  computed  from  Equations  (5.  3)  and  (5.4).  The  problem  is 
to  determine  the  extreme  pressure  levels  for  the  same  frequency  band 
at  point  3. 


flow  ^  #2  #2 


£  1 

l2 

3  2  r 


Figure  7.  Two  Pipes  in  Series. 


Equation  (3.  8)  gives 


Pg  =  ea2^2 


1  +  +  2A2  cos 

1  +  T"*  +  2  r  cos  (Ig 
1  2  2 


(5.  5) 


Unlike  the  termination  impedance,  Z ^  cannot  be  considered  constant 
over  the  frequency  band.  In  fact,  its  rate  of  change  with  respect  to 
frequency  is  great.  Thus,  for  all  but  very  narrow  bandwidths,  f  0 
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ranges  from  zero  to  1,  and  @2  will  range  from  zero  to  ir  in  the  fre¬ 
quency  band. 

The  maximum  Pg  obviously  occurs  under  the  following  condi¬ 
tions: 

1.  Pg  is  a  maximum. 

2.  cosipy  =  1  . 

'  3  .  co  s  6 2  =  -1  .  , 

4.  J"*2  is  a  maximum. 

Similarly,  minimum  Pg  will  occur  under. .the  conditions: 

1.  Pg  is  a  minimum. 

2.  cos  \p2  =  -1. 

3.  cos  02  =  !• 

4.  is  a-  maximum. 

It  can  be  shown  in  both  cases  that  conditions  1,2,  and  4  cannot 
occur  at  the  same  frequency.  Furthermore,  the  possibility  of  condi¬ 
tion  2  occurring  simultaneously  with  even  one  of  the  other  three  condi¬ 
tions  is  very  remote. 

It  was  hoped  that  even  with  the  stated  discrepancies  the  condi¬ 
tions  might  give  levels  which  would  provide  an  estimation  of  the  actual 
extremes.  Example  2  of  Chapter  VI  discounts  this  idea  by  showing 
the  levels  to  be  completely  unreasonable.  Successive  application  of 
the  conditions  to  series  systems  with  more  than  two  pipes  magnifies 
the  factor  of  unreasonability. 

In  conclusion,  it  seems  the  analytical  unpredictability  of  the 
junction  impedance  functions  and  6  relative  to  the.  function  i p  defies 
analytical  determination  of  maximum  and  minimum  allowable  pressure 
levels  for  complex  systems. 


CHAPTER  VI 

EXAMPLES 


6,  1  Introduction. 

Application  of  piping  system  analysis  methods  developed  in  Chap¬ 
ters  I-V  is  illustrated  by  three  numerical  examples.  The  examples  are 
oriented  toward  the  problem  of  specification.  The  first  example  deals 
with  open  end  impedance  calculation.,  Determination  of  maximum  and 
minimum  allowable  pressure  levels  for  a  series  piping  system  is  at¬ 
tempted  in  the  second  example.  The  third  example  illustrates  the 
method  of  complete  analysis  for  a  relatively  simple,  but  practical, 
piping  system. 

The  pound-foot-second  system  of  units  is  used  throughout.  All 
pressure  amplitudes  are  expressed  in  decibels  by  the  conversion 

db  =  20  log10  Tp“  (6.1) 

o 

where  P  denotes  pressure  pulsation  magnitude  in  psf  and  the  reference 
pressure  PQ  =  4.  18  x  10  ^  psf  (.  0002  dynes  per  cml  ).  The  fre¬ 
quency  range  of  interest  is  0  to  10,  000  cps  in  each  example. 

The  specified  termination  pressure  amplitude  in  the  second  and 
third  examples  is  taken  as  zero  decibels  over  the  complete  frequency 
range.  Results  of  the  examples  can  then  be  added  to  any  other  speci¬ 
fied  level  in  decibels  to  obtain  allowable  levels  for  the  desired  speci¬ 
fication. 
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The  liquid  considered  has  the  following  physical  properties: 
p  =  2.0  slugs/ ft 3.  ; 
p  =  2.4x10  8  lb-sec/ft2. ; 
c  =  4900  ft/ sec. ; 

K  =  4.  75  x  107  psf. 

9 

The  elastic  modulus  of  the  pipe  material  is  2  x  10  psf.  The 
constant  n  is  taken  to  be  1.8  for  all  pipes. 

6.  2  Example  1.  Open  End  Impedance. 

The  purpose  of  this  example  is  to  numerically  show  the  frequency 
characteristics  of  the  open  end  impedance  discussed  in  Section  4.  3.  The 
impedance  functions  R^(x)  and  Xj(x)  defined  by  Equations  (4.  4)  and 
(4.  5)  are  shown  graphically  in  Figure  8  as  reproduced  from  Kinsler 
and  Frey 

For  digital  computation  purposes,  approximating  polynomials  for 
0  <  x  <  8  were  calculated  from  numerical  data  in  Kinsler  and  Frey. 

These  polynomials  are 

R X(x)  =  (0.  124324)  x2  -  (5.  051  x  10"3)  x4  + 

(9.  618  x  10' 5)  x6  -  (9.  199  x  10"7)  x8  +  (6.  2) 

-Q  If) 

(3. 621  x  10  ) x  U  , 

and 

Xx(x)  =  (0.  4220302)  x  -  (2.  73694  x  10"2)  x3 + 

(7. 08794  x  10'4)  x5  -  (8. 45110  x  10'6) x7  + 

(3.  9 1 7 5 0  x  10"8)  x9  . 

Each  polynomial  has  a  standard  error  of  estimate  of  0.  0016  for  28  points. 

Equations  (6.  2)  and  (6.  3)  were  used  in  Equation  (4.  2)  to  calculate 
the  open  end  impedance  of  a  5.  157"  inside  diameter  pipe.  The  reflection 
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coefficient  of  the  open  end  was  computed  from  Equation  (3.  5)  where 
the  velocity  of  propagation  in  the  pipe' was  taken  to  be  4200  ft/ sec. 
Results  are  shown  in  Figure  9.. 

6.3  Example  2.  Maximum  and  Minimum  Analysis. 

The  series  piping  system  shown  in  FigurelO  is  used  to  illustrate 
the  methods  of  finding  allowable  maximum  and  minimum  pressure  levels 
discussed  in  Chapter  V.  Pertinent  system  data  is  provided  by  the  fig¬ 
ure.  Pipe  diameters  are  assumed  such  that  the  velocities  of  propaga¬ 
tion  are  as  given. 

The  input  impedance  was  calculated  in  the  previous  section.  For 
this  example,  the  reflection  coefficient  was  assumed  constant  over  one- 
third  octave  filter  bands  as  shown  in  Figure  llf  Using  these  constant 
values,  P2max  ancl  ^2 min  were  computed  from  Equations  (5.  3)  and 
(5.  4)  where  Pr  is  taken  as  zero  decibels.  Pressure  levels  at  the  points 
3  and  4  were  calculated  on  the  basis  of  conditions  outlined  in  section 
(5.4). 

Results  of  the  analysis  are  given  in  Figures  12a  and  12b.  As 
expected, the  allowable  levels  at  points  3  and  4  are  unreasonable. 

6.4  Example  3.  Complete  Analysis. 

Figure  13  depicts  a  suction  line  of  a  piping  system.  System  data 
is  given  on  the  figure.  The  analysis  objective  is  to  determine  the  allow¬ 
able  pressure  amplitude  spectrum  at  the  pump  so  that  the  pressure 
amplitude  spectrum  of  zero  decibels  at  the  terminus  is  not  exceeded. 

The  strainer  plate  at  the  termination  poses  a  somewhat  different 
situation  than  previously  discussed.  It  is  evident  that  adding  the  hole 
impedances  as  parallel  impedances  will  give  the  total  impedance  Z  . 
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=  4200  ft/ sec..  Example  2. 
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Center  Frequency  (1/3  Octave  Bandwidth) 

Figure  12b.  Maximum  and  Minimum  Allowable  Pressure  Levels  at  Points  3  and  A,  Example  2. 


Figure  13.  Suction  Line  to  Illustrate  Complete  Analysis  Procedure,  Example 
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If  the  hole  impedances  are  identical,  is  found  by  dividing  the  impe¬ 
dance  of.  one  hole  by  the  total  number  of  holes.  The  problem  is  then  to 
determine  the  impedance  of  each  hole.  In  this  sample,  each  hole  of  the 
strainer  plate  is  assumed  to  have  the  same  impedance  which  is  given  by 
the  method  of  open  end  impedance  determination  discussed  in  Section 
4.  3.  The  physical  situation  violates  the  infinite  baffle  condition  imposed 
on  this  method,  but  this  discrepancy  does  not  seem  to  merit  a  more  in¬ 
volved  study  at  this  time.  The  transition  in  area  from  the  strainer  plate 
to  the  first  pipe  is  neglected. 

Actual  computations  needed  for  analysis  are  described  by  the 
FORTRAN  program  shown  in  Figure  14.  Necessary  explanations  ac¬ 
company  the  program.  Relationships  coming  either  in  full  or  in  part 
from  equations  in  preceding  chapters  are  noted.  Necessary  revisions 
of  these  equations  due  to  mean  flow  in  the  direction  of  positive  lateral 
direction  have  been  made  as  outlined  in  Section  2, 4.  Erroneously,  the 
velocity  of  propagation,  c,  was  taken  as  4100  ft/sec  in  the  calculations. 

For  qualitative  results,  computations  for  freoaency  increments 
not  to  exceed  4  cps  were  found  to  suffice.  For  frequency  ranges  of  spe¬ 
cial  interest,  the  increment  of  computation  was  reduced  to  0.  5  cps. 

The  allowable  pressure  spectrum  at  the  pump  was  found  to  be 
somewhat  periodic.  A  typical  segment  of  the  spectrum  is  shown  in 
Figure  15.  The  critical  part  of  this  segment  is  in  the  frequency  range 
4500  <f<  4750  cps.  It  was  found  that  the  spectrum  contained  waveforms 
almost  identical  to  the  critical  part  of  this  segment  every  1050-1100 
cps;  the  first  such  waveform  occurring  at  approximately  500  cps. 

It  was  also  found  that  all  "local"  absolute  maximums  and  mini- 
mums  occurred  in  these  critical  parts  of  the  spectrum.  Thus,  to 
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obtain  the  maximum  and  minimum  allowable  pressure  levels  sought 
in  Chapter  V,  calculations  are  needed  only  for  these  critical  ranges  of 
frequency.  Figure  16  shows  the  maximums  and  minimum s  obtained 


in  this  manner. 


52 


DIMENSION  ZC(3),V(3},a(3),ALF(3),$INH(3),COSH(3) 


Pl-3. 1415926 

PI  »  * 

READ  1  .CONI, C0N2 

C0N1  [s  (r)  2]^rr  con?  'if [4  (ej] 

DO  7  1-1,3 

READ  l,ZCO),V(l).,EL(l),ALF(1) 

ZC(I>  =  Zc.  .  V(I)  *  Oj  ,  EI.(I)  ■  tt  ,  ALF(I)  =  ,, 

ARG-ALF(I)*EL(1) 

ARG  * 

E-EXP(ARG) 

SINH(  1  )■(£•*  1 .0/E)/2.0 

SINH(I)  =  sinh  tfjij 

C0SH( l)-(E+1.0/£)/2.Q 

COSH  (I)  =  cOBb  <ti»l 

IF  (SENSE  SWITCH.  2)  4,3 

READ  1.C.DC.CM 

C  *  f(cpsl  ,  DC  »  Af  ,  CM  =  fma<<  . 

GO  TO  5 

ACCEPT  1,C 

Accept  arbitrary  f(cpn)  from  console  typewriter 

WiC*2,0*PI 

W  =  u  =  2*f 

X2-C0N1*C**2 

X2  =  Re  |"lmpedance  at  the  seaj  »  l.  09  x  10'^  f2 

Y2-CON2*C 

Y2  »  Im  ^Impedance  at  the  seaj  =  0.663  f 

1-1 

Pipe  Number  1 

EXECUTE  PROCEDURE  10 

X1-SH*ZC( »') 

Xl.(Zc|)neJsinhT1',J 

Y1-$HJ*ZC(I) 

Y1  =  (Z0|)Im(sir.l.T,/|-) 

EXECUTE*  PROCEDURE  2 

'  J 

EXECUTE  PROCEDURE  30 

G 1 DB-20 . 0*0 . 43429440*LOG ( GAB ) 

GIDB  =  gain  in  db  across  pipe  no.  1 

XI -CH 

Xl  «  Re  ^ cosh  j 

Yl-CHJ 

Yl  =  Im  /cosh  y  .t.  1 

EXECUTE  PROCEDURE  2 

1  J 

X2«X3-SH/ZC(I) 

Y2-Y3-SHJ/ZC(l) 

Xl-G 

YI-GJ 

EXECUTE  PROCEDURE  2 

Z1-X3 

Z1  3  I}e  £  Impedance  at  junction,  pipe  no.  1  J 

Z1J-Y3 

Z1J  =  Im  ^Impedance  at  junction,  pipe  no.  1^ 

1-2 

Pipe  Number  2 

EXECUTE  PROCEDURE  10 

X 1«ZC( 1 )*CH 

Y  l-ZC( 1 )*CHJ 

X2-SH 

Y2-SHJ 

EXECUTE  PROCEDURE  2 

Z2--X3 

Z2  «=  Re  ^Impedance  at  junction,  pipe  no.  2^ 

Z2J— Y3 

Z2J  =  Im  /impedance  at  Junction,  pipe  no.  2  J 

XI-Z1*Z2-Z1J*Z2J 

1 

Yl-Z l*Z2J+Z2*ZI J 

r  “  lropedance  at  junction,  pipe  no.  3 

X2-Z1+Z2 

Y2-Z1J+Z2J 

EXECUTE  PROCEDURE  2 

Eq.(4.2) 


Eq.(*.2) 


Kq,  (3. 3V 


Eq.  (6.  1) 
Bq.(4.2) 


I2q,  (3.  1) 


Eq.(4.n 


Eq.  {3.  n 


Figure  14.  FORTRAN  Program  for  Complete  Analysis  Procedure, 
Example  3. 
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EXECUTE  PROCEDURE  10 
X1-SK*ZC(l) 

Y1»SHJ*ZC{ I ) 

X2-X3 

Y2-Y3 

EXECUTE  PROCEDURE  2 
EXECUTE  PROCEDURE  30 
G2D8-20 .0*0 .  43429448*L0G.(GAB ) 
G»G I DB+G2DB 

IF  (SENSE  SWITCH  I.)  20.21 

20  PUNCH  1 3,G1DB,G2DB,G,C 
GO  TO  22 

21  TYPE  13,G1DB,G20B,G,C 

22  IF  (SENSE  SWITCH  2)  4,23 

23  C-C-fOC 

IF  (C-CM)  5,5,40 

BEGIN  PROCEDURE  -2 

Q«X2*X2+Y2*Y2 

X3«(XI*X2+Y1*Y2)/Q 

Y3-(X?*Y1-X1*Y2)/Q 

END  PROCEDURE  2 

BEGIN  PROCEDURE  30 

G-CH-X3 

GJ-CHJ-Y3 

ARG-G*G+GJ*GJ 

GAB-SQR(ARG) 

END  PROCEDURE  30 
BEGIN  PROCEDURE  10 
B«W*EL( I )/V( I ) 

SB-SIN(B) 

CB-COS(B) 

CH-COSH( I )*CB 
.  CHJ-SINH(I)*SB 
SH-S 1 NM( I )*CB 
SHJ-COSH ( I )*SB 
END  PROCEDURE  10 

1  FORMAT  (E10.3,E10.3,E10.3,EI0.3) 
13  FORMAT  (F7.I ,F7.l.F7v1,3X  F8.1) 
END 


Pipe  Number  3 

Xl-(Zc3)  R« 

Y1  ■  (Zc3)  to  |.inh  y,l3  } 

X2  =  Re,  ^Impedance  at  junction,  pipe  no.  3 
Y2  *  Im  ^mpedance  at  junction,  pipe  no.,  3 


G2DB  *  gain  in  db  acroaa  pipe  no.  3 
G  *  system  gain  in  db 


Increment  frequency 
Teat  for  maximum  frequency 

(Xi)  *  j(Y3)  ■  :  j|ysi 


GAB  ««  J  Preaaure  gain) 


B  - 

SB  =■  sin  (Ujij) 


CB  =  coa  Ojij) 

CH  =  Re  ^cosh  y^/j^ 
CHJ  *  Im  ^cosh  yjf ^ 
SH  =  Re  {slnh  yjljJ 


SHJ  =  Im  /sinh.y. 


J‘> 


^  Input,  output  format  control 


Eq.  (3.  3) 


Eq.(3,  l) 


Eq.(6.1) 


Eq.(3.  3) 


Figure  14.  FORTRAN  Program  for  Complete  Analysis  Procedure, 
Example  3. 


Figure  15.  Typical  Segment  of  Allowable  Pressure  Spectrum  at  Pump,  Example  3. 
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